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Abstract 

We study the synchronization transition (ST) of a modified Kuramoto model on two different 
types of modular complex networks. It is found that the ST depends on the type of inter-modular 
connections. For the network with decentralized (centralized) inter-modular connections, the ST 
occurs at finite coupling constant (behaves abnormally). Such distinct features are found in the 
yeast protein interaction network and the Internet, respectively. Moreover, by applying the finite- 
size scaling analysis to an artificial network with decentralized inter-modular connections, we obtain 
the exponent associated with the order parameter of the ST to be (5 ~ 1 different from /?mf ~ 1/2 
obtained from the scale-free network with the same degree distribution but the absence of modular 
structure, corresponding to the mean field value. 

PACS numbers: 89.75.Hc, 05.45.Xt, 89.75.-k 
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Recently, there have been considerable efforts to understand complex systems in terms of 
graph, consisting of vertices and edges Q. An interesting feature emerging in such complex 
networks is a power-law degree distribution, Pd{k) ~ A;" 7 , where degree k is the number of 
edges incident upon a given vertex. Such networks are called scale-free (SF) networks |2|. 
Many SF networks in real world such as the yeast protein interaction networks (PINs), the 
Internet, and social networks contain modular structures within them While 
recent studies have focused on the identification or classification of such modules Q, little 
attention has been paid to how such modules are in—ted eaohoth.B.In^te, 
we study the synchronization transition (ST) emerging from such modular SF networks. It 
is found that the feature of the ST depends on the types of the inter-modular connections, 
different from that occurring on the SF network without modules. 

Synchronization phenomena can be found in diverse fields 0, , which have been studied 
through coupled oscillator models [11]. Here, we study the dynamics of synchronization on 
SF networks through a modified Kuramoto model 
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where 0, and u>i are the phase and the intrinsic frequency of vertex i, respectively. u>i is 
chosen from the Gaussian distribution with unit variance. The summation runs over j, 



the nearest neighbors (nn) of vertex i. Note that in contrast to other works 11 
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coupling strength is used as K/ki, where k{ is the degree of vertex i. Due to the factor 
l/k{, the dephasing effect at the vertex with large degree is reduced and the ST occurs 
at finite K c in SF networks even for 2 < 7 < 3 
absent, K c = for SF networks with 2 < 7 < 3 
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defined as M. = lim 
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. Note that when the factor 1/h is 
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141 . 1151 ] . The order parameter A4 is 



where (■ • • ) is the ensemble average over different 



re total number of vertices. M. is (1) in the fully incoherent 



configurations and iV is t 
(coherent) phase. 

To study how the profile of inter-modular connections affects the ST pattern, we first 
investigate the inter-modular connections of two real world networks, the yeast PIN and 
the Internet at autonomous system level. Although the two networks are not oscillator 
networks, we use them because they are prototypical examples of modular networks and 
contain different types of inter-modular connections. While a neuronal network is a good 
example of oscillator networks, its structure has been still veiled. As shown in Fig. Q the 
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FIG. 1: (Color online). Profile of inter-modular connections for the PIN (a) and the Internet (b). 
The concentration Cj,- = £i j£j i/LiLj, where li j is the number of edges emanating from module i 
to module j and the Li = NA, SA, AF, OC, EU, and AS stand for the North America, 

South America, Africa, Oceania, Europe, and Asia, respectively. 

modules of the PIN are interwoven diversely, while those of the Internet are connected 
mainly through only a single module, corresponding to the North America continent. The 
synchronization transition patterns occurring on those two networks appear distinct in the 
three aspects. First, as the coupling constant K increases, the order parameter increases 
very drastically for the former, while it does gradually for the latter (Fig.|2J). Quantitatively, 
the order parameter increases from A4 = 0.2 to 0.8 by the increment 5K = 1.3 and 5K = 
2.9 for the PIN and Internet, respectively. Second, the susceptibility \i defined as x — 
\J ((r — (r)) 2 ) with r = lim^oo \(1/N) X^Li ! |> exhibits a peak at the transition point 
K c , which is narrow (broad) for the former (latter) (the inset of Fig. EJ). To quantify it, we 
measure the interquartile range jl(| , which corresponds to the standard deviation often used 
for the asymmetric distribution case, to be 1.3 and 3.0 for the PIN and Internet, respectively. 
Each number means the interval of K where 50% of the data around the peak position belong 
to. Third, when K ^> K c , for the former, most modules are synchronized with almost the 
same phase, while for the latter, individual modules are synchronized independently, leading 
the overall system to be weakly synchronized. We also confirm such different types of the 
synchronization patterns through two artificial networks below. 

Two artificial networks are defined as follows: One is a modification of the community 
model (CM), proposed by Girvin and Newman |2[. There, n c modules are present from the 
beginning, which is fixed under evolution, while the total numbers of vertices and edges grow 
with time. At initial, each module contains m vertices and they are fully connected. At each 
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FIG. 2: Plot of the order parameter A4 as a function of the coupling constant K for two different 
types of complex networks, the yeast PIN(O) and the CM (□) in (a) and the Internet (O) an d the 
hierarchy model (□) in (b). Insets : Plot of the susceptibility as a function of K for each case. 

time step, a new vertex is added to the system and chooses a module at random. The newly 
added vertex attaches m edges to existing vertices selected with the probability proportional 
to its degree [2J, belonging to the chosen module. This process is repeated until the step 
when each module contains no vertices on average. Then n c sets of distinct SF networks are 
generated. After those steps, we add m e edges, connecting vertices belonging to different 
modules, which are selected randomly. The parameters, n c , no, mo and m e are adjustable. 
To compare the resulting network with the yeast PIN, we take n c = 5 and no = 625, mo = 3 
and m e = 934, making the mean degree be 6.58, comparable to that of the PIN, 6.55. While 
the CM may be too artificial or have too many parameters, it simply reflects the presence of 
modules and their diverse connections. The second one is the hierarchy model introduced by 
Ravasz and Barabasi jl^. The model contains 5 modules at the highest level, each of which 
contains 5 submodules within it in a hierarchical manner. The inter-modular interactions 
in the hierarchy model are made only through one particular module among the 5 modules 
at the same hierarchical level. The mean degree of the hierarchy model is (k) = 7.8. 

The distinct ST patterns can be observed more obviously in the size-dependent behaviors 
of M. in the CM and the hierarchy model. The system size N is controlled by changing the 
number of vertices in each module no in the CM and the number of hierarchical levels in the 
hierarchy model. As shown in Fig. El for the CM, the ST occurs at K c ~ 2.1, which is the 
crossing point of Ai for different sizes. We find numerically that the order parameter scales 
as M ~ N~ /3 / tJ -F(AN 1 ^), where A — K — K c , (3/fi ~ 0.25 and 1/u ~ 0.25. Thus (3 w 1 
for the CM network. To check if it is related to the existence of the modular structure, 
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FIG. 3: Size dependence of A4 in the CM (a) and the hierarchy model (b). The data are obtained 
for different system sizes 125(0), 625(0)> 3125(d) and 15625(A) with 1000 configurations for each. 
Inset: Plot of the order parameter in the scaling form. 




FIG. 4: The normalized phase distribution of each module (solid line) at the coupling constant 
K = 7.0 for the PIN (a), the CM (b), the Internet (c), and the hierarchy model (d). The dashed 
line represents the normalized phase distribution of the overall system. 

we destroy the modular structure in the CM model by rewiring edges. Two different edges 
are randomly selected and their one-end vertices are exchanged unless this swap procedure 
makes the network disconnected. The swapping process is applied as many as two times 
of the total number of edges. Then the resulting network loses modules but the degree 
distribution remains the same. Applying the finite-size scaling analysis to the swapped CM 
network, we obtain the mean field value /?mf ~ 1/2- Thus the nature of the ST is changed 
by the absence of modular structure. On the other hand, for the hierarchy model, we cannot 
find any crossing point in M. for different sizes. For small K, the data are well fit to the 
scaling form of M ~ N-^F(KN X ^) with p/p w 0.5 and ~ 0.17 (Fig.Eb). 
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FIG. 5: The interplay of the ST pattern in the hierarchy model as a function of p, the swapping 
ratio, for p = (O), p = 0.1 (□), p = 0.5 (A), and p = 1.0 (o). Inset: Plot of ^ max as a function of 
P- 

To understand the emergence of the different ST patterns from a microscopic respect, 
once we identify modules of each network and examine the synchronization within each 
module. For the yeast PIN, 43 modules are identified according to their functional cate- 
gories and for the Internet, six modules are done according to the continents in the world. 
Modules for the artificial networks can be naturally identified. In Fig. HJ we plot the phase 
distributions at K = 7.0 for the four networks. Each curve represents the normalized phase 
distribution of each module. For the PIN, the phase distributions of each module overlap 
each other, but the three modules, cell adhesion, mitochondrial transcription and septation. 
For the CM, most of the phase distributions overlap with each other, implying that the 
networks are globally synchronized. For the Internet and the hierarchy model, however, the 
phase distribution functions of each module do not overlap, indicating that each module 
is synchronized independently. The network is modularly synchronized. Thus, the global 
synchronization for the Internet and the hierarchy model appears weaker than that for the 
PIN and the CM for the given K. 

To see the interplay between the two ST patterns, we deform the network structure of 
the hierarchy model by swapping edges as used in the CM network. The ratio p between 
the number of swapping edges and the total number of edges is a control parameter. With 
increasing p, the inter-modular connections are decentralized. As can be seen in Fig. the 
ST pattern becomes sharper as the ratio p increases, and in particular it changes drastically 
up to p ~ 0.1, beyond which they look very similar and overlap with the ST data for the 
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FIG. 6: The ST patterns as a function of K for various mean numbers of emanating edges of the 
BA model, m B A = 1-0 (o), 1.1 (□), 1.2 (A), 1.4 (o), and 2.0 (o). Inset: Plot of 4^ as a function 
of m B A- 



PIN or the CM. 

It was recently shown that th e sy nchronizability in complex networks is related to the 
maximum value of load £ mSLX 1$, llfll ]. Load is the effective amount of data packets passing 
through a given vertex when every pair of vertices sends and receives a unit data packet 
along the shortest pathways between them j^j]. To examine how £ max is related to the 
synchronization pattern, we investigate £ max as a function of the swapping ratio p. As 
shown in the inset of Fig. £ max decreases rapidly in the range of < p < 0.1, however, it 
decreases very slowly or is almost constant for p > 0.1. This suggests that across p ~ 0.1, 
the synchronizability becomes enhanced drastically, in accordance with the behavior of the 
M. 

It would be interesting to see such a crossover through the Barabasi and Albert (BA) 
model. We control the number of edges wba emanating from a newly added node in the BA 
model. When tuba — 1, Mis almost flat with changing K. For larger wba, M. changes more 
sharply with increasing K as shown in Fig. El Cm decreases drastically up to itlba ~ 1.1 
and gradually decrease beyond that as shown in the inset of Fig. El reflecting that the 
synchronizability enhances drastically at ttl-qa ~ l l- This fact is related to the topological 
feature as follows. When ttl-qa — 1, the network is tree. As fn^A increases, the edges 
connecting different branches are formed, and its fraction is non-negligible for tji-qa ^ 11- As 
a result, the synchronization transition occurs more drastically for m-QA > 11- It was found 
that as m-QA increases, the exponent 8 associated with the load distribution, Pi{l) ~ £~ s , 
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2l|. This fact 



increases drastically from 5 m 2.0, and almost remains as 2.2 for wba > 11 
is closely related to the crossover behavior of A4. 

In summary, we have studied the two different types of the synchronization patterns 
on SF networks through a modified Kuramoto model. Such different patterns are mainly 
rooted from different types of inter-modular connection profile. When the inter-modular 
connections are decentralized (centralized), the ST occurs very drastically (gradually) and 
the synchronization arises globally (within modules), of which the behavior was analyzed 
quantitatively. We also showed that the ST point occurs at finite coupling constant for the 
former, while it is not well identified for the latter by performing finite-size scaling analysis. 
The critical exponent (3 associated with the order parameter is obtained to be (3 m 1 for 
the former case. The crossover behavior between the two cases was also studied by rewiring 
edges from the hierarchy model and by modifying the BA model. Finally, it is interesting 
to note that the global or the modular synchronization may be related to epileptic seizure. 
While the brain is electrically activated, any disruption within it may cause abnormal func- 
tioning. Then neurons in one hemisphere misfire and generate abnormal electrical activities, 
which spread to the other hemisphere through diverse pathways, leading the global synchro- 
nization |22j. To treat it, the corpus callosum, the main communication pathways between 
the two hemispheres, is sometimes cut to prevent the communications from one side to the 
other. It reveals through recent experiment 2^ that while the synchronization after the cut 
arises within one module mostly, the global synchronization also arises occasionally, which 
is accomplished via other pathways. 

This work is supported by the KOSEF Grant No. R14-2002-059-01000-0 in the ABRL 
program. 
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